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Abstract — The probability density function (PDF) and cumu- 
lative distribution function of the sum of L independent but not 
necessarily identically distributed squared 77-/1 variates, applica- 
ble to the output statistics of maximal ratio combining (MRC) 
receiver operating over 77-/7 fading channels that includes the 
Hoyt and the Nakagami-m models as special cases, is presented 
in closed-form in terms of the Fox's H function. Further analysis, 
particularly on the bit error rate via PDF-based approach, is 
also represented in closed form in terms of the extended Fox's H 
function (H). The proposed new analytical results complement 
previous results and are illustrated by extensive numerical and 
Monte Carlo simulation results. 

Index Terms — ?/-/i variates, cellular mobile radio systems, non- 
integer parameters, diversity, maximal ratio combining, binary 
modulation schemes, bit error rate, Fox's H function, Meijer's 
G function, Fox's H function, and extended Fox's H function. 



I. Introduction 

IN recent times, different diversity schemes have marked an 
important impact in the arena of wireless communication 
systems. The main reason behind this is that these different 
diversity schemes allow for multiple transmission and/or re- 
ception paths for the same signal [1]. The optimal diversity 
combining scheme is the maximal ratio combining (MRC) 
diversity scheme where all the diversity branches are processed 
to obtain the best possible signal-to-noise ratio (SNR) [1]- 
[3]. This results into extensive occurrence of the statistical 
distribution of the sum of squared envelopes of faded signals 
in several wireless communication systems [4]. 

Additionally, wireless communications are driven by a com- 
plicated phenomenon known as radio-wave propagation that 
is characterized by various effects such as fading, shadowing 
and path-loss. The statistical behavior of these effects is 
described by different models depending on the nature of 
the communication environment. Various distributions have 
several applications in wireless communication engineering 
problems and one of those that we focus on is a communi- 
cation system employing MRC diversity scheme undergoing 
77-/1' distribution i.e. the study of MRC diversity combining 
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'in this work, both the Formats are considered and hence the results 
presented herein are applicable to both the Formats. 



receiver operating over rj-fi fading channels [5], where the 
statistics of the sum of squared 77-/7 random variates (RVs) are 
required. Moreover, the performance analysis of such wireless 
communication systems usually requires complicated and te- 
dious tasks related to statistics as explained in detail in [6]. For 
instance, it is worth sharing here that the probability density 
function (PDF) and cumulative distribution function (CDF) of 
the sum of L independent but not necessarily identical (i.n.i.d.) 
Gamma RVs have been investigated in [7] using the results 
presented in [8] to address the performance of MRC diversity 
receivers over Nakagami-?7i fading channels. Similarly, in [9]- 
[11], the authors have independently presented a closed-form 
expression for the PDF of the sum of independent Gamma 
variates with arbitrary fading parameters. 

Recently, the 77-/1 fading distribution that includes the 
Nakagami-777 and the Hoyt distribution as its special cases, 
has been proposed as a more flexible model for practical 
fading radio channels [12] and free-space optical (FSO) com- 
munication systems. Additionally, the 77-/7 model has been 
found suitable for modeling small scale fading phenomenon in 
non-line-of-sight mobile communication channels [12]. Non- 
homogeneous physical characteristics of the environment have 
been taken into consideration in the mathematical modeling of 
this channel that sometimes best fits to the practical scenario 
and experimental data [13]. The 77-/7 distribution fits well to 
experimental data and can accurately approximate the sum 
of i.n.i.d. Hoyt envelopes having arbitrary fading degrees 
[14]. In particular, its tail closely follows the true statistics 
where other distributions fail to yield a good fit [15]. Beyond 
that, this distribution has been gaining interest in the field 
of performance evaluation of digital communications over 
fading channels (for example, see [4] and references therein). 
Furthermore, in [4], the authors have published results on 
the distribution of the sum of the squared envelopes of 77- 
/i faded signals but the published results are given in terms 
of nested infinite-series summations representation, dependent 
on the number of channel branches and utilizing the confluent 
Lauricella hypergeometric function [16]. It has been stated in 
[17] and other relevant research proceedings that the confluent 
form of the multivariate Lauricella hypergeometric function is 
inherently difficult to be evaluated especially when its order 
increases. Moreover, this function does not have a simple real- 
axis integral representation. 

In this work, we offer novel closed-form expressions for 
the PDF and CDF of the sum of i.n.i.d. squared 77-/7 RVs in 
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terms of easily computable Fox's fi function [18]-[20, App. 
(A. 5)]. It is noteworthy to mention that the bit error rate 
(BER) is one of the most important performance measures 
that forms the basis in designing wireless communication 
systems. Hence, we also propose new closed-form expressions 
of the BER, as a performance metric, for binary modulation 
schemes, via a PDF-based approach, of a L-branch MRC 
diversity receiver in the presence of 77-/1 multipath fading, 
in terms of the extended Fox's fi function (H) [21] 2 . This 
proves the importance and the simplicity in the employment 
of those earlier derived simple closed-form statistical PDF and 
CDF expressions. These resulting easy-to-evaluate expressions 
also give an alternative form for previously known/published 
results. It should be noted that all our newly proposed results 
are readily computable using the Mellin-Barnes theorem that 
further corroborates the generality and the usefulness of the 
analytical frameworks introduced in this paper. These have 
been checked and validated by Monte Carlo simulations. 

The remainder of the paper is organized as follows. Section 
II introduces the system and Section III gives novel closed- 
form expressions for the PDF and CDF of the sum of squared 
77-/X RVs in terms of the Fox's fi function. Next, Section IV 
utilizes these results presented in Section III to derive useful 
expressions for the BER, as a performance metric for MRC 
diversity receivers operating over i.n.i.d. 77-/1 fading channels 
or diversity paths in terms of the extended Fox's fi function 
(H). Further, Section IV also discusses the results followed 
by the summary of the paper in the last section. 

II. Channel and System Models 

A MRC based communication system with a source and 
a destination is considered with L diversity paths undergoing 
i.n.i.d. 77-/1 fading channels as follows 



The PDF of 7/ can be expressed in two formats, described 
in what follows. 

• The 77-/1 Fading Model - Format 1: In Format 1, the 
parameter < r\i < 00 denotes the scattered-wave power 
ratio in a non-homogeneous environment between the in-phase 
and quadrature components of the fading signal within each 
cluster of multipath that are assumed to be independent from 
each other and to have different powers. In such a case, 
hi = 2+ \ 1+ '" and Hi = Vl . It is noted that in Format 
1, H/h = (l-r,)/(l + ri). 

• The 77-/1 Fading Model - Format 2: In Format 2, — 1 < 77; < 1 
denotes the correlation coefficient in a non-homogeneous 
environment between the scattered-wave powers of the in- 
phase and quadrature components of each multipath cluster. 
In such a case, hi = , 1 $ and Hi = . m a . It is noted that in 
Format 2, H/h = 77. 

It can be easily seen from above that one Format may be 
obtained from other by the bilinear relation 771 = (1 — J72)/ (1 + 
772) and/or equivalently 7/2 = (1 — r7i)/(l + 771), where 771 is 
the parameter 77 for Format 1 and 772 is the parameter 77 for 
Format 2 [15]. 

As discussed earlier, the 77-/1 fading distribution comprises 
both Hoyt (/i; = 0.5) and Nakagami-m (77; — ► 0, 7/; — ^ 
00, 77; — > ±1) as its special cases. 

In a MRC combiner, received signals from all diversity 
antennas are algebraically added together and the combiner 
is connected to a suitable detector to detect the information 
signal. The output SNR of the combiner can be expressed by 
[1] 



1MRC 



L 

1=1 



u 1=1 



(3) 



where the instantaneous SNR of the Ith branch is given by 
7; = (Eb/No) \ai\ 2 , Ei, is the average energy per bit, and Nq 



Yi = aiX + n h l = l,2,...,L, 



(1) is the one sided power spectral density of the AWGN. 



where Yi is the received signal at the receiver end, X is the 
transmitted signal, on is the channel gain, and n; is the additive 
white Gaussian noise (AWGN). In an 77-/1 multipath fading 
channel, 7/ = \ai\ 2 follows squared 77-/1 distribution. Hence, 
the channel gains experience multipath fading whose statistics 
follow squared 77-/1 distribution with PDF given by 



P~fi (7) 



-cxp 



x I„ 



2/i;fe;7i 
li 

( 2mHqi 
\ li 



(2) 



where /i; > is known as fading figure representing the 
diversity order of the fading environment, T(-) denotes the 
Gamma function [22, Eq. (8.310)], /„(.) is the modified 
Bessel function of the first kind and arbitrary order v [22, 
Eq. (8.406.1)], and 77 = E [71] is the average SNR with E [.] 
denoting expectation. In 77-/1 fading model, 2/i; represents the 
number of multipath clusters [12], [13]. 

2 The extended Fox's H function (H), the Fox's H function, the Fox's H 
function, and the Meijer's G function are extensively defined in [9], [10, Table 

q. 



III. Closed-Form Statistical Characteristics for 
the Sum of Squared 77-/1 Random Variates 

This section presents the results on the statistical charac- 
teristics including the PDF and CDF of the sum of i.n.i.d. 
squared 77-/1 variates. 

Theorem 1 (PDF of the sum of squared 77-/1 RVs). Let {^i\f=i 
be a set of i.n.i.d. squared r]-fi variates with parameters fii and 
rft. Then, the closed-form PDF of the sum 

L 

y = J2^ 1 (4) 

2=1 

can be expressed in terms of Fox's H function 3 as 



Pv{y) = H 



0,2L 
2L.2L 



exp(y) 



C7« "( 2 ) 
1-1 L ' 



(5) 



where y > 0, the coefficient sets si , Hi , \ and St f , 



k € N are defined as 



(0,4i,pi),...,(0,A,i.,p fc )> 



(6) 



3 An efficient MATHEMATICA® implementation of the Fox's H function 
is available in [9], [10, Table II] (similar to [23]-[25]). 
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C(2) 



= (0,Si,Ati) 
(-l,A U fJ>i) 



and 



n (4) 



(-l.-Bl)Ml) ! 

respectively with A\ 



2m{h l -H l ) 



. , (0,B k ,p, k ) , 

• , (-1, B k , fi k ) , 
and B 



(7) 
(8) 

(9) 



2 t i l (h l +H,)- 



Proof: In order to derive the PDF of Y, we proceed as 
follows. Firstly, the MGF 



7W 7i (s) ^E[cxp(- 7i s)] 



exp(-7;s)p 7! (7)d7 (10) 



of a single 77-/1 distribution is given as [4] 

M 7l (s) = (1 + sAi)-^ 1 (1 + sBi) 



(11) 



Then, after performing some simple algebraic manipulations 
using [26, Eq. (6.1.15)], the MGF of a single r)-fi distribution 
can be rewritten as 



r«(i + svi i ) rw ( 1 + s - B i) 

T^(2 + sAi) T^{2 + sBiY 



(12) 



Since, 7 ; 's are independent, the MGF of Y is the product of 
the MGF's of the 7z 's, i.e. 



M Y (s) = UM yi (i 



(13) 



i=i 



Now, the PDF of the sum of squared r\-p, RVs, using the 
obtained MGF in (13), via inverse Laplace transform [27] can 
be expressed as 



PY (y) = C- l {M(s)} 



l 

2tt7 



M Y (s)exp(s)ds (14) 



that produces a Mellin-Barnes contour integral [28] represen- 
tation as 



Pv{y) 



2m 



r«(l + sAi) r«(l + sBi) 
c --rw(2 + S^;)r«(2 + s5;) 



exp(s)<is. 
(15) 

Hence we use this obtained result and perform some simple 
rearrangements on the terms in the Mellin-Barnes contour 
integral representation to express the closed-form PDF of the 
sum of squared r)-fj, RVs, Y, in terms of the Fox's fi function 
as given in (5). ■ 

It is worthy mentioning that the PDF of sum of squared 
77-/1 RVs has also been successfully achieved in terms of the 
confluent form of the multivariate Lauricella hypergeometric 
function in [4] that typically involves an L-fold infinite sum- 
mations. On the other hand, our alternative result presented in 
Theorem 1 involves only one single-fold integration. 

The motivation and possibility that led to the above result 
presented in Theorem 1 was the representation of the PDF of 
a single squared 77-/1 RV in terms of the Fox's fi function as 
discussed below in Corollary 1. 

Corollary 1 (PDF of a single squared 77-/1 RV). Let {7;} 
be a squared 77-/1 variate with parameters fii and rji. Then, 



the closed-form PDF of this single squared r\-p RV can be 
expressed as 



PiAli) = H 



0.2 



exp(7; 



(0 > i4j, W ),(0,-Bi,w) 

(-i,Ai,m), (-i,Bi,m) 



(16) 



where A\ and Bi are as given under Theorem 1. 



Proof: The PDF of a single squared 77-/1 RV can be 
expressed using the obtained MGF in (12), via inverse Laplace 
transform [27] that produces a Mellin-Barnes integral [28] 
representation. Hence we use this obtained result to express 
the PDF of a single squared 77-/1 RV in an alternative form, in 
terms of the Fox's fl function, as expressed in (16). ■ 
Fig. 1 presents the PDF of the output SNR obtained from the 
exact closed-form expression (5) and shows a perfect match 
between this obtained closed-form analytical result and the one 
obtained via Monte Carlo simulations for varying L's (i.e. L = 
3, 4, 5), and their respective fixed fading parameters = 1, 
/12 = 1.5, /13 = 2, /14 = 3.5, and /15 = 4.5. The value of 77 
was fixed at 1.2 representing Format 1. Similar results have 
also been obtained for Format 2. 



between Analytical and Simulation Results 




Output Signal-to-Noise Ratio y 



Fig. 1 . Comparison between PDFs obtained analytically and via Monte Carlo 
simulations for varying branches L and respective fixed fading parameters for 
these channels as fii = 1, ^2 = 1.5, fJ.3 = 2, 114 = 3.5, and /15 = 4.5. 



Theorem 2 (CDF of the sum of squared rj-fi RVs). The CDF 

ofY can be closely expressed in terms of the Fox's H function 
as 



iV(y) = l + S! 



0,2L+1 
2L+1,2L+1 



exp(y) 



s«s£\ (1,1,1) 

Bf.sW (0,1,1) 



(17) 

where the coefficient sets Sjj. , ^\ and St. 4 ^ are defined 
earlier in (6), (7), (8), and (9), respectively. 

Proof: In order to derive the CDF of Y, we proceed as 
follows. 

We integrate the PDF expressed in (16) from through 7 
and obtain the CDF for a single squared 7/-/1 RV, in terms of 
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the Fox's H function, as 



Comparison between Analytical and Simulation Results 



^(7/)=H 



exp(7;) 



(0,A l ,p l ),{0,B l ,p l ), (1,1,1) 
{-l,Ai, iH ),{-l,B hf xi), (0,1,1) 



(18) 

where A\ and Bi are as given under Theorem 1 . 

Now, performing a similar integral operation on (5), utiliz- 
ing a similar explanation as presented in the proof of the PDF 
of the sum of squared r)-/j, RVs i.e. Theorem 1 to obtain (5) 
from (13), and further making some simple modifications to 
the Mellin-Barnes integral representation to satisfy the exact 
definition of the Fox's fi function, we obtain a final closed- 
form result for the CDF of Y, in terms of Fox's fi function 
as presented in (17). Hence, in other words, the expression 
presented in (18) is a special case of (17) in Theorem 2 for 
L= 1. ■ 

Fig. 2 and Fig. 3 present the CDF and the logarithmic CDF 
respectively of the output SNR obtained from the exact closed- 
form expression (17) and show a perfect match between this 
obtained closed-form analytical result and the one obtained via 
Monte Carlo simulations for varying L's (i.e. L = 2,3,4,5), 
and their respective fixed fading parameters /ii = 1, fis = 1.5, 
/i3 = 2, fii = 3.5, and /15 = 4.5. The value of ?/ was fixed 
at 1.2 representing Format 1. Similar results have also been 
obtained for Format 2. 



Comparison between Analytical and Simulation Results 




Output Signal-to-Noise Ratio y 



Fig. 2. Comparison between CDFs obtained analytically and via Monte Carlo 
simulations for varying branches L and respective fixed fading parameters for 
these channels as fii = 1, = 1 .5, = 2, /14 = 3.5, and ^5 = 4.5. 



IV. Applications to the Performance of Diversity 
Combining Receiver Systems 

This section presents the results on the performance anal- 
ysis, in particular outage probability (OP) and BER analysis, 
of the sum of i.n.i.d. squared r\-\x variates. In here, we limit 
ourselves to binary modulation schemes. 

In MRC combining scheme, all the branches are selected at 
the output. In our case, for a L-branch MRC diversity receiver, 
the SNR y, is given by (3). 




Fig. 3. Comparison between CDFs obtained analytically and via Monte Carlo 
simulations, on log scale, for varying branches L and respective fixed fading 
parameters for these channels as /^i = 1, fj,2 = 1-5, ^3 = 2, /I4 = 3.5, and 
M5 = 4.5. 



A. Outage Probability 

When the probability that the instantaneous MRC output 
SNR falls below a given threshold y t h, we encounter a 
situation labeled as outage and it is an important feature to 
study OP of a system. Hence, another important fact worth 
stating here is that the expression derived in the Theorem 2 
also serves the purpose for the expression of OP of MRC 
diversity combining receivers based wireless communication 
system that is experiencing i.n.i.d. 77-yU fading channels or in 
other words, when the desired user is subject to rj-/j, fading, 
the probability that the SNR falls below a predetermined 
protection ratio yth can be simply expressed by replacing y 
with y t h in (17) as 



-Pout (yth) = -fV(yth)- 



(19) 



Employing similar substitutions, all the other respective ex- 
pressions of CDF can be utilized for OP such as replacing 7; 
with 7th in (18) and/or replacing y with y tn in (17). 

B. Average BER 

The most straightforward approach to obtain BER P e for 
MRC is to average the conditional BER 4 for the given 
SNR, over the PDF of the combiner output SNR i.e. P e = 

J °°P e (e\y)pY (y)dy [1]. 

Theorem 3 (BER of a L-branch MRC system operating over 
r\-[L fading channels for binary modulation schemes). The BER 
of a L-branch MRC diversity combining receiver wireless 
communication system running over r/-fi fading channels, 
valid for any binary modulation scheme including coherent 



4 The expression P e (e|j/) 



unified conditional BER ex- 



r(p,g«) 

pression for coherent and non-coherent binary modulation schemes over an 
AWGN channel [29]. r(-, ■) is the complementary incomplete gamma function 
[22, Eq. (8.350.2)]. The parameters p and q account for different modulation 
schemes. For an extensive list of modulation schemes represented by these 
parameters, one may look into [30] or refer to Table I. 
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TABLE I 

Conditional Bit Error Probability (BER) Parameters 



Modulation p q 



Coherent Binary Frequency Shift Keying (CBFSK) 0.5 0.5 

Coherent Binary Phase Shift Keying (CBPSK) 0.5 1 

Non-Coherent Binary Frequency Shift Keying (NBFSK) 1 0.5 

Differential Binary Phase Shift Keying (DBPSK) 1 1 



q = 0.5. In Monte Carlo simulations, the r\-p fading channel 
was generated based on the description given in [12]. 

We observe from Fig. 4 that the analytical results provide a 
perfect match to the MATLAB simulated results and the results 
are as expected i.e. the BER decreases as the SNR increases. 
Its important to note here that these values for the parameters 
were selected randomly to prove the validity of the obtained 
results and hence specific values based on the standards can 
be used to obtain the required results. 



binary frequency shift keying ( CBFSK), non-coherent binary 
frequency shift keying (NBFSK), coherent binary phase shift 
keying ( CBPSK), and differential binary phase shift keying 
(DBPSK), can be expressed in closed-form, in terms of the 
extended Fox's H function (H) 5 , as 

(l-<7,l,p),Ci 



Comparison between Analytical and Simulation Results 



P - <7_fr2L+l,l 
r e ~ 2 xl 2L+2,2L+2 



£2,(-q,l,P) 



where 



and 



Ci = T«,Ti 2) , (1,1,1), 



C 2 = (0,l,l),Tf\T ( , 4) 



(20) 

(21) 
(22) 



L ' L ' 

and where the coefficient sets , , andT k 4 \ k G N 
are defined as 



T 



(2) 



T 



(3) 



(2,A,Mi), 
(2,Bi,Mi), 
(1, A\,n\), 



..,(2,A k ,p k ), 
■ ■ j (2, B k ,p, k ), 
..,(l,A k ,p k ), 



and 



T[ 4) = {l,B 1 ,p l ),...,(l,B k ,p k ), 



(23) 
(24) 
(25) 

(26) 



respectively with Ai and Bi are as defined in Theorem 1. 

Proof: Utilizing BER equation P e by substituting the 
conditional BER equation P e (e\y) and (5) into it and per- 
forming some simple manipulations along with some simple 
rearrangements of T(.) function terms, we get an exact closed- 
form result of the integral valid for any binary modulation 
scheme including CBFSK, NBFSK, CBPSK, and DBPSK, in 
terms of the extended Fox's H function (H), as presented 
above in (20), Theorem 3. ■ 

C. Results and Discussion 

The numerical results for BER of MRC diversity combining 
receiver scheme with L-diversity over i.n.i.d. r\-p fading 
channels are presented in this section. 

The average SNR per bit in all the scenarios discussed is 
assumed to be equal. In addition, different digital modulation 
schemes are represented based on the values of p and q 
where p = 0.5 and q = 1 represents CBPSK, p = 1 and 
q = 1 represents DBPSK, CBFSK is represented by p = 0.5 
and q = 0.5, and NBFSK is represented by p = 1 and 

5 The extended Fox's H function (H) was first introduced in [21] and has 
a MATHEMATICA® implementation given in [9], [10, Table III]. 




Signal-to-Noise Ratio (dB| 



Fig. 4. Average BER of different binary modulation schemes over i.n.i.d. r\- 
[i fading channels with L = 5-branch MRC and fading parameters for these 
channels as fj,i = 1, fj,2 = 1.5, fi3 = 2, (14 = 3.5, and ^5 = 4.5. 



Furthermore, it can be seen from Fig. 4 that, as expected, 
CBPSK outperforms the other modulation schemes and the 
coherent binary modulation schemes outperform their respec- 
tive non-coherent and/or differential binary modulation scheme 
i.e. CBPSK outperforms DBPSK and CBFSK outperforms 
NBFSK. Additionally, PSK in general performs better than 
FSK, as expected. Similar results for any other values of p's 
and 77's can be observed for the exact closed-form BER for 
L-diversity i.n.i.d. rj-fj, channels presented in this work. 

V. Concluding Remarks 

We derived novel closed-form expressions for the PDF 
and the CDF of the sum of i.n.i.d. squared r]-p RVs. An 
interesting finding is that these expressions can be written in 
terms of special functions, specifically the Fox's H function 
distributions. Based on these statistical formulas obtained, and 
following the exact PDF-based approach, we analyzed the 
performance of a MRC diversity combining receiver based 
wireless communication system operating over i.n.i.d. r\-p 
fading channels and important performance metrics such as OP 
and BER were expressed in closed form and hence this serves 
as the key feature along side the novel statistical derivations 
of PDF and CDF. 

An exact closed-form expression for the BER performance 
of different binary modulations with L-branch MRC scheme 
over i.n.i.d. rj-p fading channels was derived. The analytical 
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calculations were done utilizing a general class of special func- 
tions including the Fox's H function [31], the Fox's H function, 
and the extended Fox's H function (H). In addition, this work 
presents numerical examples to illustrate the mathematical 
formulations developed in this work and to show the effect of 
the fading severity and unbalance on the system performance. 
Our results complement previously published results that are 
in the form of infinite sums or nested sums of the fading 
parameter. 
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